Coulomb blockade in graphene quantum dots 
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Abstract 

We study the conductance spectrum of graphene quantum dots, both single and multiple cases. 
The single electron tunneling phenomenon is investigated and the periodicity, amplitude and line 
shape of the Coulomb blockade oscillations at low temperatures are obtained. Further, we discuss 
the transport behavior when multiple dots are assembled in array and find a phase transition of 
conductance spectra from individual Coulomb blockade to collective Coulomb blockade. 
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Introduction. After first isolation from bulk graphite, graphene has attracted intense 
experimental and theoretical attention due to its unusual electronic spectrum and hence 
exotic properties [l|, 0]. Moreover, its zero nuclear spin and weak spin-orbit interaction 
make it one of the most promising platforms for solid state quantum information processing. 
Design and fabrication of graphene quantum dot to confine Dirac fermion should be a critical 
step towards realizing such promise. There are alternative theoretical schemes proposed to 
deal with graphene's gapless spectrum and form quantum dot Recently, single 

or double graphene quantum dots by etching nanoribbon into a charge island have been 
demonstrated 0, 0, 0, [I3] • Although there are some initial studies 11, 12], the transport 



properties of graphene quantum dot are still remain to be explored both experimentally and 
theoretically. 

In this paper, we study conductance of graphene quantum dots in both single dot and 
multiple dots regime, respectively. In previous theoretical studies, the conductance reso- 
nances are determined by the level of the (quasi-)bound states of the noninteracting elec- 
trons 0, [ll|, 12]. In the experiment the shape of the resonances will be governed by the 
electron interaction via the Coulomb blockade (CB) effect . We use the Constant In- 

teraction (CI) model to account for the intradot Coulomb interaction and perform an exact 
numerical calculation of the low bias conductance within quantum transport formalism 13 . 
The obtained CB peaks provides a simple way to describe the rich phenomena of transport 
physics in a variety of graphene quantum dots. In particular, we find the interdot tunneling 
would lead to a phase transition from individual CB to collective CB in the conductance 
spectra through multiple graphene quantum dots. 

Graphene quantum dot. We begin by considering the graphene quantum dot on electro- 
statically confined graphene nanoribbon with semiconducting armchair boundary |4]. The 
electron wave function in graphene is usually described as a four-component (iso)spinor 

^ = {^^\ \ \ ^ b Close to the Dirac point, i.e. on low energy scales, a weakly 
doped (or undoped) graphene can be described by two identical and decoupled Dirac Hamil- 
tonian 

H K = v F o ■ p + U(x,y), (1) 
H K > = v F a* -p+U(x,y), (2) 

where the Fermi velocity v F ~ 10 6 m/s, the momentum operator p = —ih~(d x ,d y ) T , Pauli 
matrices cr = [a xi a y ), a* = (a x , — o- y ), and U(x,y) = U(y) stands for the potential of 
applied electric field along y direction. We solve the equation H^fr = respectively in 
left barrier, middle dot and right barrier to calculate the bound state energy within the 
conduction band if Vb arr ier > Vg ate - According to the semiconducting armchair boundary 
conditions, the quantized transverse momentum is given as 

(n+ 1/3W n , „ 
q n = K - ^fAn = 0,±l,--- (3) 

where W is the width of the ribbon. Additionally, after matching the wave function at 
y — and y — L (L is the length of dot), we can obtain discrete energy levels because of 
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FIG. 1: Calculated energy levels in the graphene quantum dot for different quantized transverse 
momentum (a) qo, (b) q-\. We consider a relatively long dot qoL = 5 as in Ref. 4j. The horizontal 
axis stands for AV = Vbarrier ~ Vgate an d the vertical axis represents the electron bound energy, 
both in a unit of hvpqo- 



the electrostatical confinement in the longitudinal direction. The results are systematically 
shown in Fig. 1. The bound state energy levels only exist when E satisfies 



E - eVg a te > hv F q n > \E - eV b 



barrier 



(4) 



Coulomb blockade behavior. In the following analysis, we consider the five lowest quan- 
tized energy levels in the graphene quantum dot at potential barrier AV^ = 1 in a unit of 
hvpqo, of which three belong to go, and two belong to Here we set the nanoribbon 



width W = 45 nm according to Ref. [16|, then Hvpq = 0.015 eV. First, we adopt the 
Constant Interaction (CI) model 13J], which assumes that the Coulomb interaction between 
the electrons is independent of the number N of electrons in the dot and can be described 
by a constant capacitance C, and estimate the charging energy e 2 /C ~ 1 in a unit of hvpq , 
which is constant for a given graphene quantum dot geometry and in the order of exper- 



imental measurements [7|, [l6j. In this model, the additional energy is e 2 /C + AE, where 



AE is the energy difference between two consecutive states. The CB peaks will appear if 
the gate charging can supplement the additional energy for next electron. Then we use the 
method described in Ref. 13j to discuss in detail the single electron tunneling phenomenon 
at low temperature in the weak tunneling regime (hT <C &bT AE). The linear response 
conductance is 



G 



2k B T 



J2^n f'(E N + U(N) - U(N - 1) - E F ) 



(5) 



N 



where f(x) 



l+e k B J 



is the Fermi-Dirac distribution function, E^ is the energy of the 



top filled single electron state for a A^ electron dot, U(N) = (Ne 2 )/2C — Ne<p ext in which 
(fiext = 0o + ctV ga t e) and E F is the chemical potential of the lead. The width of localized 
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FIG. 2: Coulomb blockade oscillations for a graphene quantum dot at different temperatures. 
From top to bottom, k B T = 0.02,0.06,0.10,0.14,0.18,0.22,0.26,0.30, in a unit of hv F qo- Spin 
freedom has been considered. The gate voltage is in a unit of eV and all the conductance peak 
heights are normalized by the first peak at fcgT = 0.02. 

energy T in the above equation is determined by tunneling though the classically forbidden 



region 



14 



T 



exp — 



L+d 



k 



dy 



(6) 



where k! = i^/q 2 — ((E — eVbarrier) / 'hv p) 2 is the vanishing wave vector in the barrier region 
and d is the width of each barrier (setting q d = 1). Notice that the tunneling rate will be 
suppressed at a large k'. 

Fig. 2 shows the calculated conductance of graphene quantum dot as a function of gate 
voltage at different temperatures, in which we take consideration of spin freedom. We find 
that the derived function of Fermi-Dirac distribution performs as a Delta function and it 
induces a conductance peak at position E N + U(N) — U(N — 1) = E F . The peak height 
maximum G max , according to Eq. (5), is given by G max = gjj-fF, which decreases linearly 
with increasing temperature in the quantum CB regime. Also, the iVth peak probes the 
specific excitation spectrum around En, the quantum CB regime therefore usually shows 
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15 . We would also like to mention that in recent 



randomly varying peak heights 0, [HI, 
experiments [7|, la, |9] graphene dots have most likely rough edges which can induce electron 
scattering at the atomically sharp edges. Nevertheless we can find the lack of rough edges 
effect in the model described here does not change qualitatively the main features of the 
CB. Furthermore, the line shape of conductance can be obtained as 

G , i, 5 



cosh-'(— — ), 5 = E N + U(N) - U(N - V 



Ei 



(7) 



Gmax QikfiT 

which is consistent with the recent experimental observations in single electron transistor 
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on graphene nanoribbon 16 . 

Graphene quantum dot array. In the following section, we investigate the conductance 
spectrum multiple coupled graphene quantum dot systems realized by assembling the above 
single dot separated from each other by a metallic gate. The type of coupling between 
quantum dots determines the character of the electronic states and the nature of transport 
through the quantum dot array. Here we use a Hubbard model to account for the effects of 
nearest interdot tunneling which is much smaller than the Coulomb energy of each dot, and 
neglect the long-ranged Coulomb interactions which are screened by the metallic gate 

We consider a doubly spin-degenerate state in each graphene quantum dot and to study 
the tunneling phenomenon in coupled dots. Due to the weak spin-orbit coupling in graphene, 
the spin degeneracy cannot be easily broken if without applying magnetic field. The tun- 
neling strength in this case is determined by the interdot distance exponentially If the 
distance between the nearest dots is quite large compared with the dot size, the tunneling 
strength approaches zero, and the individual dot energy levels are conserved. As the interdot 
distance is decreased, the tunneling between nearest dots cannot be ignored, and the energy 
states are split into subbands, while spin degeneracy is still remained. As the discrete energy 
levels including the intradot Coulomb repulsion have been calculated above, it is convenient 
to get the tunneling coupled energy levels via direct diagonalization 
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Fig. 3 shows the calculated conductance of two, three and four dots, from which, we find 
a remarkably phase transition for tunneling coupled graphene quantum dot array [17|. For 
weak interdot tunneling, i.e., when D/L is large (D is the distance between nearest dots 
and L is the length of a single dot), the Coulomb blockade of individual dots is maintained. 
The number of peaks is equal to the number of energy levels within a single graphene dot 
without interacting with surrounding dots. Each peak represents the addition of electrons 
to the array, one to each graphene dot at the same time. For the moment, transport 
through double graphene dots in this regime has been reported 10]. For intermediate D/L, 
the Coulomb blockade of individual dots is destroyed and a collective Coulomb blockade 
phenomenon appears. The original two energy states are split into two subbands and each 
capacitance peak is split into several peaks equal to the dot number. For stronger interdot 
tunneling strength, subbands will cross each other and the Coulomb blockade is destroyed 
altogether. All these theoretical predictions can be verified in the near future experiments. 

Conclusion. In summary, we investigate the quantum dot behavior of electrostatically 
confined graphene nanoribbon. We study the single electron tunneling phenomenon at low 
temperatures of this system. The periodicity, amplitude and line shape of the CB oscillations 
are discussed in detail. Also, we have presented calculations of conductance spectra for 
multiple coupled graphene dots and find a phase transition from individual CB to collective 
CB as the coupling strength increases. The results presented here are important to provide 
necessary information for future experimental work. 
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(c) 

FIG. 3: Plots showing the evolution of the Coulomb blockade oscillations as a function of the 
interdot tunneling strength D/L. D is the distance between two nearest graphene dots and L is 
the length of individual dot. Here we demonstrate the conductance spectrum of multiple coupled 
graphene dots for the case of (a) two, (b) three, (c) four. The gate voltage is in a unit of eV and 
the conductance is normalized by the highest peak. 
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